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Abstract 

The spreading of vector-transmitted diseases, like dengue, is a relevant problem of public health. 
Different techniques are used to understand its complex dynamics. In this paper we use a stochastic 
Markovian dynamics approach to describe the spreading of dengue and the threshold of the dis- 
ease. The coexistence space is composed by two structures representing the human and mosquito 
populations. The human population follows a susceptible-infected-recovered (SIR) type dynamics 
and the mosquito population follows a susceptible-infected-susceptible (SIS) type dynamics. The 
human infection is caused by infected mosquitoes and vice- versa so that the SIS and SIR dynamics 
are interconnected. After setting up the master equation associated to the stochastic dynamics 
we develop a truncation scheme to solve the evolution equations from which we get the epidemic 
curve and the threshold of the disease as well as the reproductive ratio. The results are compared 
to the actual outbreaks of the disease. 

PACS numbers: 87.10.Mn, 87.10.Hk, 05.70.Fh, 05.70.Ln 



I. INTRODUCTION 



Dengue is a vector-borne infectious disease with very complex dynamics, whose spreading 
is a relevant problem of public health. The disease is transmitted to human mainly by 
the mosquito Aedes Aegypti. Many factors are determinant for the transmission of dengue 
in urban centers such as the climatic conditions for the vector proliferation, the human 
concentration and mobility. Although many efforts are put in the development of a vaccine 
against the four types of virus, until now the only available strategy to reduce the spreading 
of the disease pQ is the control of the vector population. Therefore, it is very important 
to analyze the effect of vector control avoiding the occurrence of dengue epidemics. In this 
context, the inter-host modeling of dengue dynamics and control may be a very useful tool 
for helping its understanding and the establishment of vector control strategies. 

Different techniques and approaches are used to model the dynamics of transmitted dis- 
eases [21 [3] such as deterministic differential equations [4], stochastic dynamics [5H7], cellular 
automata [8], and complex networks [§]. Concerning vector transmitted diseases including 
dengue modeling, there are also different schemes and approaches [21 E] based on deter- 
ministic models of differential equations [T0HT2] , probabilistic cellular automata [13HT6] and 
complex networks [T71 IT8] . However, as far as we know, it is the first time that a description 
in terms of a master equation defined on a lattice is used on the investigation of a vector- 
borne infectious disease although it was used by some of us to investigate epidemic models 
of direct transmitted diseases [5] . 

As in other works [T0TTT21 [15], the present approach is also motivated by actual data of 
dengue epidemics [19] in particular, by two outbreaks of dengue occurred in Salvador, Bahia, 
Brazil in 1995 (without vector control) and 2002 (with vector control). Those data had also 
motivated a previous analysis based on the basic reproductive ratio [TO] . 

The first dengue model was proposed by Newton and Reiter [20] in 1992 assuming a 
susceptible, exposed, infected and recovered (SEIR) structure for humans and susceptible, 
exposed and infected (SEI) structure for mosquitoes due to the fact that the mosquitoes die 
before being removed. This framework has been followed by other continuous and discrete 
dengue models [10HT21 [Til EE]- Here we consider a simpler model, illustrated in figure [TJ 
assuming a susceptible, infected and recovered (SIR) structure for humans and susceptible, 
infected and susceptible (SIS) structure for mosquitoes. The infection of humans is due to 
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FIG. 1: Illustration of the reactions. Left panel: SIR structure. A susceptible human (Hg) 
becomes infected (Hi) through a catalytic reaction mediated by infected mosquitoes (Mi). An 
infected human becomes recovered (Hr) spontaneously and remains permanently in this state. 
Right panel: SIS structure. A susceptible mosquito (Ms) becomes infected through a catalytic 
reaction mediated by infected humans. An infected mosquito spontaneously becomes susceptible. 

mosquitoes and the infection of mosquitoes is due to humans. In other words the infection 
reactions, S— >1, on both structures are catalytic and not autocatalytic as happens to the 
original SIR and SIS models [5H7j. The other reactions, I— >R in the SIR structure and I— >S 
in the SIS structure, are spontaneous reactions. The mosquito structure has been simplified 
by suppressing the death and the birth of mosquitoes which amounts to saying that a dead 
mosquito is immediately replaced by a new-born (susceptible) mosquito. We are, therefore, 
assuming that the number of mosquitoes remains constant throughout the outbreak of the 
epidemics. 

The features of our model that differ from that of Newton and Reiter [20] are as follows. 
Firstly, we do not distinguish between susceptible and exposed states both for humans and 
mosquitoes. Secondly, deaths of humans are not considered since the human life time is 
much larger than the period of the disease. As to the deaths of the mosquitoes, they are 
implicit in the model in the following sense. The reaction / — >■ S for the mosquitoes is to be 
interpreted as the death of an infected mosquito and the simultaneous birth of a susceptible 
mosquito. The major difference however rests on the use of a stochastic lattice model, which 
takes into account the spatial distribution of humans and mosquitoes. 

This paper is organized as follows. In section [IT] we introduce the model and derive from 
the master equation the evolution equation for the densities. We also define the quantities 



that characterize the spreading of the epidemics. In section III we develop the simple mean- 
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field analysis. In section IV we introduce the evolution equation for the two-site probabilities 
and perform a pair mean-field analysis. Also in this section we compare the results coming 
from the pair mean-field with actual data coming from the outbreak of the dengue disease 
occurred in Salvador in 1995 and 2002. Finally, concluding remarks are presented in the last 
section. 



II. MODEL AND EVOLUTION EQUATIONS 

The modeling of disease spreading that we consider here corresponds to a continuous 
time stochastic Markovian process defined on a lattice where the sites are occupied by 
human individuals or by mosquitoes. In order to properly describe the human and mosquito 
populations we consider two sublattices of the whole lattice, one for each population. The 
sublattices, named H and M, are interpenetrating in such a way that the nearest neighbor 
sites of a site of one sublattice belong to the other sublattice. Each site of the sublattice M 
can be in one of two states, either occupied by a susceptible mosquito Ms or by an infected 
mosquito Mj. Each site of the sublattice H can be in one of three states, occupied by a 
susceptible human H$, occupied by an infected human Hj or occupied by a recovered human 
Hr. The system evolves in time according to the following stochastic dynamics. 

At each time step a site of the lattice is chosen at random. If the site belongs to sublattice 
H then the following update is used, (a) If the chosen site is occupied by a susceptible 
human then it becomes infected with rate a times the fraction of infected mosquitoes in its 
neighborhood. If the chosen site is occupied by an infected human then it becomes recovered 
spontaneously with rate c. Once recovered the individual remains permanently in this state 
which means that, if the chosen site is occupied by a recovered human, it remains unchanged, 
(b) Next, we consider the case in which the chosen site belongs to the M sublattice. If it 
is occupied by a susceptible mosquito then it becomes infected with a rate b times the 
fraction of infected humans in its neighborhood. If the chosen site is occupied by an infected 
mosquito it spontaneously becomes susceptible with rate e. Our simple model is therefore 
described by four parameters a, b, c and e. In the applications we will further simplify by 
setting b = a. 

For convenience we introduce a stochastic variable r]i associated to each site % of the lattice 
that takes the values 0, 1, 2, 3, or 4 according to whether the site i is occupied, respectively, 
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by a susceptible mosquito (Ms), an infected mosquito (Ms), a susceptible human (H$), an 
infected human (Hj) or a recovered human (Hr). The time evolution of the probability 
distribution P(rj) of configuration i] = {rji} is governed by the master equation 

jP(v) = J2{MAv)P(Av) - Mv)p(ji)Y (i) 

i 

where Wi(rj) is the transition rate from rj to rj = Ait] and A% is the operator that changes 
rji — >■ as follows: 0— >1, 1— >-0, 2—7-3, 3— >4; and is the inverse of A{. The transition 
rate Wi(rj) is defined according to the rules stated above. The time evolution of an average 
(fiv)) = E?7 f(v)P(v) is obtained from the master equation and is given by 

j t (f(v)) = E({f(^v)-f(v)H(v)}. (2) 

Instead of using the full probability distribution P(rj) we consider an equivalent descrip- 
tion in terms of the various marginal probability distribution, obtained from P(rj). The 
time evolution of the several marginal probability distributions are obtained from the mas- 
ter equation ([I]) and comprises a hierarchic set of coupled equations which is equivalent to 
the master equation. This approach is convenient because it allows to obtain a solution of 
the set of equations by a truncation scheme to be explained shortly. In what follows we 
assume translation invariance and isotropy. 

Let us denote by the P(r)i) the marginal one-site probability that represents the prob- 
ability that a site % is in state rji and by P(rji,rjj) the marginal two-site probability that 
represents the probability that site i is in state rji and a neighboring site j is in state rjj. 
Other marginal probabilities are denoted in an analogous way. The evolution equation of 
the one-site probability P(l), which represents the density of infected mosquitoes, can be 
obtained from equation ([21) if we recall that P(l) = (Sfa, 1)) and is given by 

|p(l) = 6P(03)-eP(l). (3) 

where we used the relation P(03) = (S(rji, 0)5(rjj, 3)). The equation for P(0), the density of 
susceptible mosquitoes, can be obtained from equation pi) by using the property 



P(0) + P(1) = 1. (4) 

The evolution equations for P(2) and P(3), the densities of susceptible and infected 
humans, respectively, are obtained similarly from equation ^ and are given by 

|P(2) = -aP(12), (5) 



d 

db 



P(3) = fl P(12) -cP(3). 



(6) 



The equation for P(4), the density of recovered humans, can be obtained from equations ^ 
and (Jsl) by taking into account the property 



P(2) + P(3) + P(4) = 1. 



(7) 



To characterize the threshold of the epidemic it is convenient to write equations (J3]) for 
the density of infected mosquitoes and ^ for the density of infected individuals in the form 

|p(l)=6P(0|3)P(3)-eP(l), (8) 

|p(3) = aP(2|l)P(l)-cP(3), (9) 

where P(2|l) = P(12)/P(l) is the conditional probability of occurrence of a susceptible 
individual given an infected neighboring mosquito and -P(0|3) = P(03)/P(3) is the condi- 
tional probability of occurrence of a susceptible mosquitoes given an infected neighboring 
individual. Using the simplified notation x = P(l) and z = -P(3) the set of equations Q 
and (|9~1) can be written as 

1 -e 6P(0|3)\ 1 
oP(2|l) -c j 

where P(0|3) and P(2|l) are the cross transmission probabilities. At the early stages of 
the epidemic the cross transmission probabilities P(0|3) and P(2|l) can be considered to be 
constant (independent of time) and the set of equations (JIJ) and ^ becomes a linear set 
of equations. A fundamental quantity that characterizes the spreading of the disease is the 
so-called reproductive ratio Po which is defined here as follows 

ab 



dx/dt^ 



dz/dt j 



V 




(10) 



Pn 



ce 



-P(2|1)P(0|3). 



(11) 



The threshold of epidemic is determined by the largest eigenvalue of the matrix (10), which 
is related to the reproductive ratio by 



Pn 



or by 



-{-(e + c) + A /(e-c) 2 + 4ecJ2 }. 



(12) 



(13) 



According to the linear analysis the threshold of epidemic occurs when the largest eigen- 



value vanishes which happens, according to equation (12), when Rq = 1. Moreover, when 
A < there is no transmission of disease. According to equation (12) this happens when 
Rq < 1. The spreading of the disease occurs when A > 0, that is, when Rq > 1. The 
reproductive ratio characterizes not only the threshold of the disease but also its strength. 

III. SIMPLE MEAN-FIELD APPROXIMATION 

The three equations for the time evolution of -P(l), P(2) and -P(3) are not a set of 
closed equations because they depend on the two-site probabilities -P(12) and P(03). To 
get closed equations we firstly use a scheme, called simple mean-field (SMF) approximation, 
which consists in writing a two-site probability as the product of one-site probabilities, 
that is, P(f]i,r]j) = P(r] i )P(r]j). Using the abbreviations -P(l) = x, which is the density of 
infected mosquitoes, -P(2) = y, which is the density of susceptible individuals, and -P(3) = z, 
which is the infected individuals. This approximation gives -P(12) = P(1)P(2) = xy and 
P(03) = P(0)P(3) = [1-P(1)]P(3) = {l-x)z so that equations Q, @ and Q are reduced 
to the following forms 

^ = b{l — x)z — ex, (14) 

ft = ~ aXy > (15) 
~ = axy — cz. (16) 

To describe the spreading of the disease we consider an initial condition such that all 
individuals are susceptible so that y = 1 and z = 0. Moreover we consider a very small 
density of infected mosquitoes, that is, x = e where e << 1. Then we look for solutions for 
x, y and z in the limit t — > oo. In this limit the density of infected mosquitoes and infected 
individuals vanish, x — > and z — > 0. If the density of recovered individuals p = 1 — y — z 
approaches a nonzero value than the disease spreads. On the other hand p approaches zero 
(when e — > 0) the disease does not spread. This means that the stationary solution x = 0, 
z = 0, y = 1 (p = 0) is the solution corresponding to the non-spreading regime. To obtain 
the stability of this solution we linearize the above equations around this solution to obtain 

^ = bz — ex, (17) 

LLC 
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dp 
~dt 



ax, 



dz 
dt 



ax — cz. 



(18) 
(19) 



A linear analysis of stability amounts to determine the eigenvalues of the matrix composed 
by the coefficients of the right hand side of these equations, given by 

( -e b \ 

a • (20) 
y a — c j 

One eigenvalue is zero and the others are the roots of 



A 2 + (e + c)A + ec - ab = 0, 



that is, 



A± = ~{-(e + c) ± v /(e-c) 2 + 4a6}. 



(21) 



(22) 



The solution is stable as long as A + < 0, that is when ec > ab. The threshold of the 
spreading occurs when 

nh 

1. (23) 



ab 

ec 



For a = b, the threshold line is described by 

ec 



(24) 



as shown in the phase diagram e/c versus a/ col figure [2^,. 

The threshold obtained by the above linear analysis can equivalently be obtained from the 



condition Rq = 1, that is, when the reproductive ratio Rq given by equation (11) equals the 



unit. In the present case of simple mean-field approach P(2|l) = P(2) and P(0|3) = P(0) 
The initial condition gives, P(2) = 1 and P(0) = 1 so that 

ab 



Ro 



ce 



(25) 



which coincides with the condition (23) when Rq = 1. It is worth mentioning that Rq 



depends on the parameters a, 6, c and e only through the ratios e/c, b/c and a/c. Figure [2 
shows the lines at constant Po m the plane e/c versus a/c for the case b = a as well as the 
threshold of epidemic characterized by Rq — 1. 
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FIG. 2: Phase diagram in the plane e/c versus a/c showing the regions of spreading (S) and non- 
spreading (NS) of the disease for the (a) SMF approximation and (b) PMF approximation. The 
continuous curve is the threshold of spreading of disease. When e — >• 0, the value of a/c approaches 
zero for the SMF and the value a/c = 0.5 for the PMF. The curves inside the S region are lines of 
constant reproductive ration Rq, from top to bottom: Rq =1.2, 1.5, 2, 4, and 10. The threshold of 
spreading of disease (continuous line) corresponds to i?o = 1- 

To obtain the time dependent properties, including the epidemic curve we have solved 
the set of equations numerically starting from the initial condition x = e, y = 1 and z = 0, 
where e = 10~ 6 . An epidemic is characterized by the epidemic curve denned as the number 
of cases occurring in unit time, for instance, in a day or in a weak, plotted as a function 
of time. In other terms it is the number of susceptible individuals that are being infected 
per unit time. Instead of number of individual we may define the epidemic curve as the 
density of individuals that are being infected per unit time ( = —dy/dt. Figure [3^, shows 
an example of the epidemic curve obtained from the SMF approximation together with the 
density of infected mosquitoes. 

For sufficiently large values of time the density of recovered individuals p approaches its 
maximum value. Notice that the final density of recovered individuals p depends on a, c 
and e only through the ratios a/c and e/c and is the integral of the epidemic curve. Indeed, 
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if we integrate ( = —dy/dt, we get 

rod 

/ Cdt = p, (26) 
Jo 

where we used the initial condition y(0) = 1 and the result that the final density of infected 
individuals vanish, that is, z(oo) = 0, and the constraint p + y + z = 1. The final density 
of recovered individual, which is a measure of the size of the epidemics, may be understood 
in this approach as the order parameter in the sense that it vanishes in the non-spreading 
regime and is nonzero in the regime where the disease spreads. 

Near the threshold of epidemics the density of infected mosquitoes is much smaller than 



unit and we may neglect x in the first term on the right hand side of equation (14). The 
evolution equation for x becomes then 



dx 
dt 



bz — ex, 



(27) 



which together with equations (15) and (16) allows us to determine explicitly the size of the 



epidemics p. This is possible because these equations implies a conservation law obtained 
as follows. We start by defining the quantity <fi = cx + bz. From equations (27) and (16) its 
evolution equation is given by 

drh 

-cex + abxy. (28) 



dt 



The ratio of equations (28) and (15) gives 

d<fi 
dy 

which can be integrated to give 



ce 



ay 



-b 



(29) 



ce 



cx + bz = — lay + 6(1 — y) 



(30) 



which is the desired conservation law. The constant of integration was obtained by remem- 
bering that at t = 0, x = e — > 0, y = 1 and z — 0. 
When t — >oo, x — 0, z = and 1 — y = p so that 



ln(l - p) + Rq P = 



(31) 



where we used the relation Rq = ab/ce. This equation determines the size of the epidemics 
p and may be written in the form 



l-p 



-Rqp 



(32) 
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This equation is the same equation obtained by Kendal [21] for the model introduced by 
Kermack and McKendrick [22J to describe the directed transmitted epidemics [21 E] and also 
obtained by means of the simple mean- field approach to the SIR model on a lattice [7]. Near 
the threshold it is given by p = 2(P — 1). 



IV. PAIR MEAN-FIELD APPROXIMATION 



A. Evolution equations 

Next we set up equations for the two site probabilities. The evolution equation for the 
probability P(03) of a site being occupied by a susceptible mosquito and a neighboring site 
by an infected individual is given by 

^P(03) = eP(13) - (c + rb)P(03) + (1 - r)aP(021) + (1 - r)6P(303), (33) 

where r = 1/k, and k is the coordination number of the lattice. The notation P^, r^-, %) 
stands for the joint three-site probability where i and k are nearest neighbor sites of j. The 
evolution equation for the probability P(12) of a site being occupied by an infected mosquito 
and a neighboring site by a susceptible individual is given by 

-^P(12) = -(e + ra)P(12) + (1 - r)6P(302) + (1 - r)aP(121). (34) 

UjL 

The equations for the other two-site probabilities, P(02), P(04), P(13), and P(14) are not 
necessary because of the following relations 

P(02) + P(03) + P(04) = P(0), (35) 

P(12) + P(13) + P(14) = P(1), (36) 

P(02) + P(12) = P(2), (37) 

P(03) + P(13) = P(3), (38) 

P(04) + P(14) = P(4). (39) 

Let us consider now the equations for the time evolution of P(l), P(2), P(3), given by 
equations ^ and (JS]), and the equations for time evolution equations of P(03) and 



P(12), given by equations (33) and (34). These five equations are not closed equations 
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because (33) and (34) depend on three-site probabilities. To get a set of closed equations 



we use now a truncation scheme at the level of two-site probabilities (23]. To this end we 
start by writing the three-site probability as P(r]i,r]j,r)k) = P(r]i,r)k\r]j)P(r)j) where i and k 
are distinct neighbor sites of site j. Next we use the following approximation for the con- 
ditional probability P{r)i,r)k\r)j) = P(Vi\Vj)P(Vk\Vj)- Since p (Vi\Vj) = p (VhVj)/ p (Vj) and 
P(j]k\r]j) = P(r]k,r]j)/P(T]j) we get the following approximation for the three-site probability 



P(j]i,r]j,T]k) = P(r}i,r}j)P(rjj,r)k)/P(rjj). These results are used in equations (33) and (34) to 
get a set of closed equations in the variables P(l) = x, P(2) = y, P(3) = z, P(03) = u and 
P (12) = v. With this approximation the model is described by the set of five equations 

~- = bu — ex, (40) 
dt v 1 

I - («) 

— = av — cz } (42) 
du (v — v\v 

— = e (z - u) - (c + rb)u + (1 - r)a— — + (1 - r)b , (43) 

at y 1 — x 

civ u(ij — 

— = -(e + ra)v + (1 -r)b— + (1 - r)a — . (44) 

dt 1 — x y 

To obtain the threshold of the spreading of the disease we analyze the stability of the 

solution x = 0, y = 1, z = 1, u = and v = 0, which characterizes the state where the 

disease does not spread. To get the stability of this equation we linearize the evolution 

around this solution to obtain 

— = —ex + bu, (45) 

(JjL 

!=«,, (46) 

— = —cz + av, (47) 
at 

— = ez — (e + c + rb)u + (1 — r)av, (48) 
aL 

dv 

— = (1 — r)bu — (e + ra)v, (49) 

where p = 1 — y — z. A linear analysis of stability amounts to calculate the eigenvalues of 
the matrix composed by the linear coefficients of the right hand side of the above equations, 
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FIG. 3: Dengue outbreaks according to SMF approximation (a) and PMF approximation (b). In 
both cases the parameters are a/c = 10 and e/c = 49 and the initial condition is x = 10" 6 . The 
figures shows the number of individuals that are being infected per unit time ( = —dy/dt, or the 
density of new cases, and the density of infected mosquitoes x as functions of time. The area below 
the epidemic curve ( equals the density of recovered individuals po- For the SMF, po = 0.80, and 
for the PMF, p = 0.50. 

given by 
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Two eigenvalues are Ai = and A 2 = — e. The others are the roots of 

- (A + c)(A + e + c + rb)(X + e + ra) + ea(l - r)b + (A + c)(l - rfab = 0. (51) 
The line of stability is obtained by setting in this equation A = 0, that is, 

- c(e + c + r6)(e + ra) + (1 - r)abe + (1 - r) 2 a6c = 0. (52) 
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As before we consider b = a and write the equation that describes the threshold line as 

(1 — r)a 2 e + (1 — 2r)a 2 c = e 2 c + ec 2 + rac 2 + 2raec. (53) 

When the coordination number k — 4 (r = 1/4) the threshold line is described by 

3a 2 e + 2a 2 c = 4e 2 c + 4ec 2 + ac 2 + 2aec. (54) 

In figure [2]d we show the line described by this equation which represents the phase transition 
between the spreading to non-spreading regime. Below the transition line the non-spreading 
solution becomes unstable giving rise to the spreading solution. 

The threshold of epidemic obtained by the above linear analysis can equivalently be 
obtained in terms of the reproductive ratio. In the present case of pair mean-field approach 
P(2|l) = P(12)/P(l) = v/x and P(0|3) = P(03)/P(3) = u/z so that 

Ro = (55) 

ecxz 

Again, we may say that Ro depends on the parameters a, b, c and e only through the ratios 



e/c, b/c and a/c. Indeed, if we substitute A given by equation (12) into equation (51 ) we get 
an equation that gives Rq in an implicit form. It is easy to see that this equation contains 
the parameters a, b, c and e only through the ratios e/c, b/c and a/c. Figure [2Jd shows the 
lines at constant Ro i n the plane e/c versus a/c for the case b = a. We also show the disease 
threshold line, given by Ro = 1, that separates the non-spreading and the spreading regimes. 
Comparing the pair and simple mean field approximation we see that the inactive region of 
the phase diagram is larger for the pair approximation than for the simple approximation 
as can be seen in figure [2j This aspect can also be observed in figure [3] which shows the 
epidemic curve for both approximations using the same set of parameters. The simple mean- 
field predicts a larger number of recovered individuals, and consequently a larger number 
of individuals that have been infected, as compared to the pair mean-field. As illustrated 
in figure |3j for the same set of parameters, the final density of recovered individuals are 
Po = 0.80 for the SMF and p = 0.50 for the PMF. 

We have solved the above set of equations using the initial condition x = 10~ 6 , y = 1, 
z = 0, u = and v = 10~ 6 . From the numerical solution we have obtained the epidemic 
curve, that is the density of individuals that are being infected per unit time, ( = —dy/dt. 
Figure [3] shows an example of the epidemic curve obtained in the PMF approximation 
together with the density of infected mosquitoes. 
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It is worth mentioning that in the early stages of the spreading of the epidemics the 
quantities x, z, u and v increase exponentially according to 



x = x e xt , z = z e xt , u = u e xt , v = v e xt , (56) 



where A is the largest eigenvalue of the matrix (50). From this behavior we see that the 
ratios v/x = P(2|l) and u/z = P(0|3) are independent of time at the early stage of the 
spreading of disease, and so is the reproductive ratio Ro. 

An important aspect of our pair approach concerns the relation of the epidemic curve with 
the conditional probabilities P(2|l) and P(0|3) which are treated exactly in this approach. 
Remember that P(2|l) is the conditional probability of the occurrence of a susceptible 
individual in the presence of a infected mosquito and P(0|3) is the conditional probability 
of the occurrence of a susceptible mosquito in the presence of a infected individual. 



B. Actual epidemics data 

The theoretical analysis performed in the previous section was motivated by actual data 
of two dengue outbreaks, as in a previous work by one of us pi)]. They occurred in Salvador, 
Bahia, Brazil in 1995 and 2002: although in first outbreak no systematic vector population 
control has been applied by the local health authorities, in the second outbreak a more 
systematic mosquito control was in course due to the establishment of Brazilian Dengue 
Control Program. As it was discussed in [TO] the transmission force for the second outbreak 
is more intense than the first outbreak. 

The dengue time series of the outbreaks are shown in figure |4| It is easy to observe 
that the 2002 outbreak (prevalence of DEN3) increases and decreases faster than the 1995 
outbreak (only DEN2). It is clear that faster increasing of 2002 outbreak is associated with 
higher transmission force acquired for DEN3. But here is no consensus whether the faster 
reduction is more related to transmission force acquired for DEN3 that quickly increased the 
herd immunity of human population pool, or to vector control actions applied in Salvador 
in 2002, or both factors [TO]. 

Both time series can be fitted by means of the simple or pair mean-field approximations 
of our simple model, assuming parameter values for which the active phase is set up. Here 
we fitted the data from the two outbreaks by means of the pair mean-field approximation, as 
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FIG. 4: Pair mean- field (PMF) solutions fitted to the dengue outbreaks in Salvador, Brazil, in 
1995 (DEN2) and 2002 (DEN3 is prevalent). For both cases we used the same set of parameters, 
a = b = 2 day -1 , e = 9.8 day -1 , c = 0.2 day -1 . These parameters amounts to a reproductive rate 
#0 = 1-38. 

shown in figure |4j For both outbreaks we used the following values for the three parameters 



of the model: a = b = 2 day -1 , e = 9.8 day -1 and c = 0.2 day -1 . Using equation (55), 
these values give a reproductive rate equals R = 1.38. The total number of individuals, 
which can be understood as the carrying capacity, was considered to be 3.2 x 10 4 for the 
1995 outbreak and 6 x 10 4 for the 2002 outbreak. 



V. CONCLUDING REMARKS 



In this work we have applied stochastic dynamics to a dengue bipartite lattice model to 
analyze the transition between epidemic and non-epidemic states in terms of the probabil- 
ity of human-mosquito and mosquito-human transmission and vector control parameters. 
Analyzing the model by means of simple and pair mean-field approximations for coordina- 
tion number k, we have shown that a phase transition takes place between epidemic and 
non-epidemic states for k ^ 2. The epidemic curve resulting from the PMF approach were 
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used to fit the the patterns of time series of actual dengue outbreaks that have occurred in 
Salvador in 1995 and 2002. 

We have presented a precise definition of the reproductive rate Rq which is appropriate for 
systems described by a stochastic dynamics, that characterizes the spreading of the disease, 
and we have relate it to the largest eigenvalue of the matrix associated to the evolution 
equations. This definition can be generalized to other types of disease transmission and seems 
to be promising in the analysis of epidemics. According to our definition, the reproductive 
rate is directed related to the conditional probability of the occurrence of a susceptible human 
(mosquito) given the presence in the neighborhood of an infected mosquito (human). At the 
early stages of the epidemic, these conditional probabilities are simply equal to the unity 
in the SMF approach but are nontrivial in the PMF, what makes this approach a richer 
description when compared to the SMF. Another quantity that characterizes the epidemic 
is p, the quantity that measures the size of the epidemic and vanishes in the nonspreading 
regime. It was also determined by means of the SMF and PMF approach. 

In future work we intend to perform Monte Carlo simulations of the bipartite lattice 
model in order to analyze the spatial-temporal patterns. They may then be compared with 
the spatial-temporal patterns of 1995 and 2002 outbreaks which have different features from 
each other. In another work we will consider a model that takes into account the possibility 
of the vaccine [21] since there is a perspective of having the vaccine against dengue available 
within the next five years. 
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